Perturbative expansion of the QCD Adler function improved by 
renormalization-group summation and analytic continuation in the Borel plane 
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We examine the large-order behaviour of a recently proposed renormalization-group-improved 

expansion of the Adler function in perturbative QCD, which sums in an analytically closed form 

the leading logarithms accessible from renormalization-group invariance. The expansion is first 

written as an effective series in powers of the one-loop coupling, and its leading singularities in 

the Borel plane are shown to be identical to those of the standard "contour-improved" expansion. 

Applying the technique of conformal mappings for the analytic continuation in the Borel plane, we 

define a class of improved expansions, which implement both the renormalization-group invariance 

and the knowledge about the large-order behaviour of the series. Detailed numerical studies of 

specific models for the Adler function indicate that the new expansions have remarkable convergence 

properties up to high orders. Using these expansions for the determination of the strong coupling 

from the the hadronic width of the r lepton we obtain, with a conservative estimate of the uncertainty 

due to the nonperturbative corrections, a s (M^) = 0.31891[J' ( q^5^ , which translates to a s {M%) = 

n 1 18A+ - 0021 
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I. INTRODUCTION 



The determination of the QCD coupling constant to 
increasing precision is one of the most important goals of 
the Standard Model of particle physics (for a review see 
[l|). The non-strange hadronic decays of the r lepton are 
an important source of information on this quantity and 
have been exploited now for a couple of decades. The 
recent calculation of the Adler function to four loops in 
massless QCD Q renewed the interest in the extraction 
of the strong coupling a s at the scale of the r mass from 
the treatment of these processes 0-[T3|- In this context, 
several modifications of the perturbative expansion of the 
relevant observables have been proposed. The main am- 
biguities affecting perturbation theory are related to the 
implementation of renormalization-group invariance and 
to the large-order behaviour of the series. The differ- 
ences between the specific ways of accounting for these 
properties are the main source of theoretical error on the 
extraction of a s (AI^). 

In a recent work, Ref. (HI, we applied to the analysis 
of the hadronic r decay width a method of improving the 
perturbative expansions in QCD by summing the leading 
logarithms accessible from renormalization-group invari- 
ance, proposed in [lfl [l6{ and developed in llllial. The 
properties of the new expansion, which has been referred 
to as "improved fixed-order perturbation theory" , in the 
complex energy plane were investigated and were com- 
pared with those of the "contour-improved perturbation 



theory" (CIPT), and the standard "fixed-order pertur- 
bation theory" (FOPT). The new expansion has the ad- 
vantage of being written in an analytically closed form, 
while CIPT, the alternative approach of implementing 
renormalization-group (RG) invariance, requires the nu- 
merical solution of the renormalization-group equation 
for the strong coupling. 

It is known that the perturbative expansions of QCD 
correlators are divergent in many physically interesting 
situations, with the coefficients growing as n\ at large 
orders n (l9M22|. Alternatively, the divergent character 
of the series is inferred from the fact that the expanded 
correlators, like the Adler function, are singular at the 
origin of the coupling constant plane [23| . These prob- 
lems are present also in perturbative QED (24[, whose 
phenomenological success is explained by the fact that 
the fine structure constant is numerically very small. By 
contrast, for a relatively large coupling like a s (M^) in 
QCD the consequences are nontrivial. 

Special mathematical techniques for divergent series 
are available, like Borel summation, which under certain 
conditions recovers the expanded function from its in- 
creasing expansion coefficients. For QCD the problem 
was studied for many years, and it is known that the con- 
ditions of Borel summability are not satisfied 0, H^, H3 
(see also the reviews [Ifi}). In particular, the Borel trans- 
form of the Adler function has singularities in the Borel 
plane, known as ultraviolet (UV) and infrared (IR) renor- 
malons, the latter producing an ambiguity in the recon- 
struction of the function. However, if one adopts a cer- 
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tain prescription (e.g., the principal value prescription), 
it is possible to exploit the available knowledge on the 
large-order behaviour of the coefficients for defining a 
new expansion, in which the divergent pattern is consid- 
erably tamed. Such an approach was proposed in [27i — 
l29l | and developed recently in 0, 0, [Tl| , using techniques 
of series acceleration based on conformal mappings and 
"singularity softening" (these techniques were also ap- 
plied by other authors, for a discussion and references 
see [HI). We recall that the method of conformal map- 
ping was introduced and applied in particle physics in 
[30( for extending the convergence region of an expan- 
sion beyond the circle of convergence and for increasing 
the convergence rate at points lying inside the circle. As 
discussed in Ref. [ll|, the method is not applicable to 
the (formal) perturbative series in powers of a s because 
the expanded correlators are singular at the point of ex- 
pansion, a s = 0, but can be applied in the Borcl plane. 
It leads to a modified perturbative expansion in terms 
of a new set of functions, which have the advantage of 
resembling the expanded correlator (the Adler function), 
in particular by sharing its known singularities in the 
coupling and the Borel complex planes. 

In Refs. [f| ID, [ll[ the method was applied to the 
two standard versions of p erturbation theory, CIPT and 
FOPT. As argued in [TjJ, the new expansions are par- 
ticularly suitable in the contour-improved version, since 
they make simultaneously the RG summation and the 
Borel large-order summation of the Adler function. De- 
tailed numerical studies [(| [ll| established the good con- 
vergence properties of the latter expansions for several 
exact models which simulate the known properties of the 
Adler function. 

In this work, we consider the large-order properties of 
the expansion discussed in [l4| , which we shall henceforth 
refer to as "renormalization-group-summed" (RGS) ex- 
pansion. We investigate the properties of this scheme in 
the Borel plane and, using the techniques discussed in 
0, B El: we define a new class of expansions, which si- 
multaneously implement the renormalization-group and 
the large-order summation by the analytic continuation 
in the Borel plane. We shall refer to these as "Borel and 
renormalization-group-summed" (BRGS) expansions. 

The plan of the paper is as follows: We briefly review 
in Sec. [II] the perturbative expansion of the Adlcr func- 
tion and its connection to the hadronic decay width of r 
lepton. In Sec. IIIII we review the derivation of the RGS 
expansion following [l4| , and show further that it can be 
expressed as an expansion in powers of the one-loop solu- 
tion of the RG equation for the coupling, which is needed 
for rendering this expansion suitable for convergence ac- 
celeration. In Sec. IIVI wc discuss the properties of the 
expansion in the Borel plane, and show that it has the 
same dominant singularities as the CI expansion. In Sec. 
IVl we define a set of new Borel and RG- improved expan- 
sions, by using the technique of singularity softening and 
conformal mappings of the Borel plane [TjJ. In Sec. IVII 
we investigate the properties of the new expansions in the 



complex energy plane and illustrate their convergence for 
the physical observable relevant for the hadronic width of 
the r lepton, usinga class of models for the Adler func- 
tion considered in [5|, 0, [Tl|, [3l|, [H| . In Sec. I VIII we report 
a new determination of a s (M^) based on the new BRGS 
expansions and in Sec. IVIIII we summarize our results 
and conclusions. 



II. ADLER FUNCTION AND THE HADRONIC 
r-DECAY WIDTH 

The Adler function plays a crucial role in the determi- 
nation of a s (M^) from hadronic r decays. The method 
is discussed in the seminal paper [33j and is reviewed in 
several recent articles 0, H, 0,11 . 

The inclusive character of the total r hadronic width 
makes possible an accurate calculation of the ratio 



r[r — > z^ T hadrons] 
T\t~ —> v T e~Vp\ 



(1) 



Of interest is the Cabbibo allowed component R T y/A 
proceeding cither through a vector or an axial vector cur- 
rent, which can be expressed theoretically in the form: 
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(2) 

where N c = 3 is the number of quark colours, Sew and 
S' EW are electroweak corrections, 5^ is the dominant per- 
turbative QCD correction, and 8*z) denote quark mass 
and higher ZJ-dimensional operator corrections (conden- 
sate contributions) arising in the operator product expan- 
sion (OPE). The decay width is suitable for the precise 
extraction of the strong coupling, since the (less-known) 
higher terms in the OPE bring a very small contribution 
to ([2]). Therefore, a fairly accurate phenomcnological de- 
termination of the QCD perturbative part S^°> is possible 

The theoretical calculation of 5^ is based on unitar- 
ity, which implies that the inclusive hadronic decay rate 
can be written as a weighted integral along the timelike 
axis of the spectral function of a polarization function. 
As shown in [33| . the analytic properties of the polar- 
ization function and the Cauchy theorem allow one to 
write equivalently this quantity as an integral along a 
contour in the complex s-plane (chosen for convenience 
to be the circle \s\ = M%). After an integration by parts, 
the quantity of interest is expressed as the contour 
integral: 

where the reduced Adler function D(s) = D^ 1+0 \s) — 1 
is obtained by subtracting the dominant term from 
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the logarithmic derivative of the polarization function, 
Z)( 1+0 '(,s) = — ,sdn^ 1+0 ' l (s)/ds, where the superscript de- 
notes the spin (33|. The perturbative expansion of D(s) 
reads 0] 

oo n 
n=l k=l 

(4) 

where a s (n 2 ) = a s (/i 2 )/7r is the strong coupling at the 
renormalization scale /x. The leading coefficients c rii i are 
obtained from Feynman diagrams, the known coefficients 
c rlj i calculated to four loops in the MS-renormalization 
scheme being (see Q and references therein): 

ci,i = 1, c 2 ,i = 1.640, c 3 ,i = 6.371, c 4 ,i = 49.076. (5) 

Several estimates for the next coefficient 05,1 were made 
recently H S, & El . 

The remaining coefficients c n ^ for k > 1 are deter- 
mined from renormalization-group invariance: the func- 
tion .Dp er t, calculated in a fixed renormalization scheme, 
is scale independent and therefore satisfies the equation 



axis. The C1PT [13, [13 is defined by the choice /i 2 
when ^ reduces to 

oo 

-Dcipt(s) = c„,i(a s (-s))" , 



— s. 



(11) 



where the running coupling a s (— s) is determined by solv- 
ing the renormalization-group equation numerically 
in an iterative way along the circle, starting with the in- 
put value a s (M 2 ) at s = —M 2 . This expansion avoids the 
appearance of large logarithms along the circle \s\ = M%. 

The expansions pop and (fTTjl coincide formally as long 
as all the terms in the series are retained. In fact, the 
coefficients c„.i are known to increase as n\, so the series 
are divergent. We shall turn to this property in Sec. IIVI 
If the series are truncated at some order N , the expan- 
sions pop and (fTTj) differ by terms of order a^ +1 , this 
being the main theoretical error in the the determination 
of a s (M 2 ) from the measured r hadronic width. 



III. RENORMALIZATION-GROUP 
SUMMATION 



-Dpert(s) 



= 0, 



which can be written equivalcntly as 

dD„ 



a I \ d-Dpert 



'pert 



31n(— s/yit 2 ) 



= 0, 



where 



(0) 



(7) 



(8) 



fc=0 



is the (3 function governing the scale dependence of the 
coupling. From ([7]) one can express c„.fc for k > 1 in 
terms of c„.i and the coefficients ft of the perturbation 
expansion ([5]). The known ft coefficients, calculated to 
four loops in MS scheme, are (see [U, HH for the calcu- 
lation of ft and earlier references): 

ft = 9/4, ft =4, ft = 10.0599, ft = 47.228. (9) 

In the "fixed-order perturbation theory" calculation of 
8(°\ the choice fi 2 = M 2 is adopted, when the expansion 
(g]) reads: 



5 FO pt(s) = £(a s (M 2 ))™[c nil (10) 

n 

+ £fcc„, fe (ln(- S /M2)) fc - 1 ]. 



fc=2 



As remarked in [361 ] . due to the large imaginary part of 
the logarithm ln(— s/M 2 ) along the circle \s\ = M 2 , the 
series (fTTJj) is badly behaved, especially near the timclikc 



As proposed in [17], Qi| , the expansion (01 can be writ- 
ten in the RGS form 



D 



RGS 



( S ) = Y,(as(v 2 )TD n (y), 



(12) 



where the functions D n (y), depending on a single variable 
y = l + /3 a s (M 2 )ln(-5/M 2 ), (13) 



are defined as 



D n {y) = }^(k-n+ l)c k ,k-n+i 



k—n 



y-i 

ft 



k — n 



(14) 



As seen from the definition, the function Dx sums the 
leading logarithms in the series Q , the function D2 sums 
the next-to-leading logarith ms, and so on. The attrac- 
tive feature pointed out in HH, is that these func- 
tions can be obtained in a closed analytical form. The 
derivation is based on renormalization-group invariance: 
after inserting the expansion p~2]) into the condition ([6]), a 
straightforward calculation leads to the following system 
of differential equations for D n (y), for n > 1: 



ft 



dy 



n-l 



t=0 



(15) 

with the initial conditions D n (l) = c„.i which follow 
from (fT4"|) . 

The solution of the system p"5]) can be found iteratively 
in an analytical form. The expressions of D n (y) for n < 
5, written in terms of the coefficients Ck,x with k < n and 
ft with < j < n — 1, are: 



Di(y) 



cx.i 
V 



(16) 
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D 2 (u) = \ (c 2 ,i +01,1^2,1) , <2 2 ,i = -^-Iny, (17) 



D 3,{y) = r (c 3 ,i + c 2 ,i d 3 ,2 + Cl,l d 3 ,i) ! ( 18 ) 



d 3 ,2 
^3,1 



A) 



- lny 



(19) 



-^(l-y + lny-ln^ + ^l-y), 



D 4 (u) = ^-(c 4 ,i + c 3i i d 4 , 3 + c 2 ,i c?4,2 + ci,i d 4: i) (20) 



C?4, 3 
^4,2 

^4,1 



-3-5- In y 

Po 

-2^ (-1 + y) + ^ (-2 - 2 lny + 3 hi 2 y + 2y) 



in 



ft 
pi 



(21) 



P1P2 



^3 (-5 In 2 y + 21n 3 y + 41ny(-l + y) + (-1 + yf) - ^ (31ny + y - 2ylny - y 



Wo 



D 5 (y) = -r(c 5j i + c 4 ,i d 5 ,4 + c 3; i rf 5 . 3 + c 2 ,i d 5 2 + c x 1 d 5jl ), 



(22) 



d$A 
d5,3 
d5,2 
d5,l 



-4^ my, 
Po 

Pi 



PI 



-2\n 2 y + \ny-y+l)+^{y-l) 

Po 



% (41n 3 y - 71n 2 y + 6(y - 1) lny + (y - l) 2 ) + 2^|i(-y 2 - 3ylny + y + 41ny) 

Po Po 



6/?o 4 

Po 3 
P1P3 



4 (6 In 4 y - 26 In 3 y + 9(2y - 1) In 2 y + 6 (y 2 - 5y + 4) In y + (y - l) 2 (2y + 7)) 



(3(y - 2) In 2 y + (2y 2 - 5y + 3) lny + (y - l) 2 (y + 3)) 

PI 



lh 



+ TT^r (4y 3 -3y 2 + 61ny (y 2 - 2) - l) + £L(y - l) 2 (y + 5) - g-(y 3 - 1) 



6po 2 



3/3 2 



Pa, 2 



A) 



r 



(23) 



The expressions of D n (y) for 6 < n < 10, which depend 
also on the coefficients c n ,i for 6 < n < 10 and ft- for 

5 < 7 < 9, ar e g iven in a somewhat different form 1 in the can be written as 
Appendix of [14j . In the numerical applications presented 
in Sec lVIl we shall use the expressions of D n (y) up to n = 
18, which can be obtained easily by solving the system * J n (y) 

dHJ) with a MATHEMATICA program. 



Using Eqs. (|TB)l - (j2"3"j) and the expressions of higher 
D n (y) derived analytically, we note that these functions 



n-l 

Cn,l + Cj,irfnj (y) 



(24) 



1 For simplicity, in H3 we presented the expressions obtained by 
inserting the known numerical values of /3j for j < 3 from (|9j, 
and setting fij = for j > 4. 



where the coefheients d n j are functions of y, which de- 
pend also on the coefheients j3j of the /3-function. One 
can check that d n j(y) vanish identically for y = 1 or if 
ft = for j>l. 

By inserting ([24)) into (|12|l . we note that the denomi- 
nators y™ can be combined in each term with the powers 
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of a s (ii 2 ), so that (fT2"]) can be written as 

oo 



n=l 



n-1 
J = l 



where 



a s (-s) = 



l + P a s {^)H-s/n 2 ) 



(25) 



(26) 



is the solution of the RG equation © to one loop at the 
scale — s, written in terms of a s (/x 2 ) = a s (/i 2 )/ir. The 
terms c Mi i in the series (|25p yield the all-order summa- 
tion of the one-loop coupling (|26[) . while the remaining 
sums yield the corrections accounting for the higher order 
terms in the expansion of the /3-function. 



IV. BOREL TRANSFORM 

In this section we discuss the properties of the expan- 
sion (f2~5"j) in the Borel plane. We start by recalling the 
standard definition [|| of the Borel transform B(u) of the 
expansion (flTI) : 



B(u) 



Eu 
Cn+1,1 _ 

n=0 



(27) 



The original function Dcipt(s) is recovered from B{u) 
by a Laplace-Borel integral. Actually, the function B{u) 
has singularities on the positive axis of the u plane, so 
the Laplace-Borel i nteg ral requires a regularization. Fol- 
lowing Refs. [H, 0, lllj we shall use the principal value 
(PV) prescription. Note that as argued in [UHil, the 
PV prescription preserves the reality of the correlators in 
the s-plane, and is therefore more consistent than other 
prescriptions with the analyticity properties imposed by 
causality and unitarity. Thus, we write 



5 C ipt(s) = / cxp ( T J B ( u )du- (28) 

Po J \Poa s {-s)J 
o 

Similarly, we can define the Borel transform Bpo( u is) 
of the FOPT expansion (flQ)) . The structure of the co- 
efficients of this expansion implies that we can write 
Bfo(u, s) as: 

oo „ ri+1 / s k — 1 



B FO (u,s) = B(u) + J2-^: -y J2 kCn + 1 ' k ln 

n=0 P n - k=2 



M r) 



(29) 

The function -Dfopt(s) is obtained from its Borel trans- 
form by 



-Dfopt(s) = -5-PV / exp (- ) B FO (u,s)du. 

Po J \Poa s {M^)J 
o 

(30) 



We introduce now the Borel transform B RG s(u, y) of the 
expansion (|25p . which can be written as 

oo n n 

B RGS (u,y) = B(u) + Y^ ■g^—i^2 c jA d n+i,j(y)- (31) 

n=0 ^° ' j=\ 

The function -Drgs( s ) is recovered by the similar 
Laplace-Borel integral 



D RGS (s) = -^PV 
Po 



| B RGS (u,y)du 7 
(32) 

written in terms of the one-loop coupling (|26[) . 

As we already mentioned, the function B(u) defined 
in (|27[) has singularities on the real axis in the w-planc, 
namely along the rays u > 2 and u < — 1 [l9|, [25[ . More- 
over, the nature of the dominant singularities can be de- 
scribed exactly: they are branch points, near which B(u) 
behaves, respectively, as 



B(u) ~(l + u) 



-71 



B(u)~(l-u/2)-^, (33) 



where the exponents 71 and 72, calculated using 
renormalization-group invariance, have known positive 
values 0, [H, EH : 



7i 



1.21, 



72 = 2.58 . 



(34) 



From (f2"9"| and (J2J) it follows that these singularities 
are present also in the Borel transforms B F o(u, s) and 
Br G s{u, y). In principle, these functions might have also 
other singularities, due to the additional infinite sums 
appearing in (f2"9"|) and (pTLj) , respectively. However, as 
we shall argue below, the dominant singularities of these 
functions, i.e. the singularities closest to the origin u = 0, 
are those at u = — 1 and u = 2 contained in B{u). 

We first present some evidence which results from the 
inspection of the next-to-leading terms in the expression 
(|3"Tj) of B RG s(u, y). Thus, from the expressions (fl"6j) - (|23"l) 
we note that 



d n ,n-i(y) = -(n - 1) my. 

Po 



(35) 



By inserting this expression in (|3ip . we obtain 
by a straightforward calculation the contribution to 
B RG s{u, y) of the term with j = n: 

B RG s(u,y)\j= n = —^uB(u) lny. (36) 



1% 



Similarly, we note that 



+ (n l)(n 2) g lnV (3?) 



where 



2 P'o 2 



^) = |ln y+ ^^(l- y ). (38) 
Po Po 
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Then we obtain by a straightforward calculation the con- 
tribution of the term with j = n — 1 in (|31[) as 

BnGs(u,y)\j=n-i = -^Jr / u'B(u')du' 
Po Jo 



Pi , 



In B(u), 



(39) 



where the integral is defined along a contour from the ori- 
gin to the point u, which does not reach the singularities 
of B(u). 

The next coefficients in the second term of (j3"Tj) exhibit 
a similar pattern: d„.„-z for 1 < I < n — 1 contains 
polynomials of the variables (1 — y) and In?/, of degree 
(I — 1) and I, respectively, with coefficients depending on 
Pj, n and I. For instance, the term proportional to ln'y 
in d n>n -i has the expression 

d n , n ^t~t±lY[ {n - k Al n ly, (40) 

L k=i p o 
bringing a contribution to ()3ip of the form 



B 



RGS 



(u,y)|j=„+l_; 



tlLA 



■yu l B(u) \n l y. (41) 



This reproduces and the second term in (j3"9")l for 
I = 1 and / = 2, respectively. 

The other terms appearing in d 7ltn -i may contribute 
also with integrals of B(w) multiplied by powers of u, as in 
P9p . Thus, in general, the second term in the expression 
(f3l"|) of -Bbrg( u ; y) is expected to contain either B(u), 
multiplied by factors which vanish in the limit y — ^ 1, or 
integrals of B(u), in which the singularities have the same 
positions as in (|33|) . but are milder. Therefore, the dom- 
inant singularities of the Borel transform (1311) concide 
with the dominant singularities of the Borel transform 
B{u) defined in (|2T|). 

One may invoke also the general argument that Mueller 
[l9j used for concluding that the dominant singularities 
of the Borel transform Bpo(u, s) defined in ([2Uf are the 
same as those of B(u). The crucial observation is that 
the positions of the dominant singularities of the Borel 
transform are determined from the behaviour of the cor- 
relators in the limit of a small coupling. 2 Since in this 
limit the running coupling a s (—s) entering (|28|) . the fixed 
scale coupling a s (Mj?) entering (f3T)|). and the one-loop 
coupling d s (—s) entering (|3"2"j) are close to each other, it 
follows that the positions of the dominant singularities in 
the it-plane of the corresponding Borel transforms, B(u), 
Bfq(u, s) and Brgs(w, y), must be the same. 



2 The connection between the behaviour of the Borel-summed cor- 
relators in the a s -plane and the position and nature of the dom- 
inant singularities of the Borel transform in the u plane is given 
explicitly in the case of a finite number of renormalons in l42l . 



V. ANALYTIC CONTINUATION IN THE 
BOREL PLANE AND NEW PERTURB ATIVE 
RGS EXPANSIONS 

As discussed above, the RGS expansion (|2"5j) is diver- 
gent, the coefficients c„.i increasing as n! at large n. In 
fact, as we shall show in the next section, the divergence 
is quite bad for expanded functions supposed to resem- 
ble the physical Adler function. A procedure to tame 
this divergent behaviour is therefore mandatory. In this 
section we shall improve the large-order behaviour of the 
RGS expansion by applying a method based on analytic 
continuation in the Borel plane, applied to the standard 
CIPT and FOPT in Refs. HEl. 

The starting remark is that the Taylor expansion (|27[) 
of B(u) is convergent only inside the disk |u| < 1, lim- 
ited by the nearest singularity at u = — 1. The region 
of convergence can be enlarged if the series in powers 
of u is replaced by a series in powers of another vari- 
able. As shown in 0, [Tl], [27} , the "optimal" variable ac- 
cording to the definition proposed in (30j is the function 
w = w(u) that conformally maps the assumed holomor- 
phy domain of B(u), i.e. the whole u plane cut along 
u > 2 and u < — 1, onto the unit disk |iy| < 1 in the 
w complex plane. The expansion of B(u) in powers of 
w is convergent in the whole complex u-plane except for 
the cuts. Moreover, this optimal mapping provides the 
fastest large-order convergence rate, compared to other 
variables that conformally map onto the unit disk only 
parts of the u plane 13011 . The detailed proof of these 
statements is given in [11]. 

An additional improvement, discussed in detail in |11| . 
is obtained by exploiting the known behaviour at the 
first singularities, presented in (|3"3")) and (|34p . The main 
idea of the procedure, denoted as singularity softening 
[43]] . is to multiply B{u) with a suitable factor S(u), such 
that in the product S(u)B(u) the dominant singularities 
are compensated or are replaced by milder singularities. 
Moreover, after compensating the leading singularities, 
one can expand the resulting function in powers of vari- 
ables that take into account only the higher, i.e more 
distant, renormalons. 

As shown in the previous section, the Borel trans- 
form Brqs(u, y) of the RGS expansion (|25|) has the same 
dominant singularities as B(u). Therefore we can apply 
the techniques of imp roving the convergence mentioned 
above. Following [ll[, we consider the functions 



wim{u) 



y/1 + ll/l - y/l - Ufl 



a/1 + U/I + y/1 - Upl 



(42) 



where I , m are positive integers satisfying / > 1 and m > 
2. The function wi m {u) maps the it-plane cut along u < 
—I and u > m onto the disk \wi m \ < 1 in the plane 
wim = wi m (u). The "optimal" mapping defined above is 
w(u) = wi2(u), for which the entire holomorphy domain 
of the Borel transform is mapped onto the interior of the 
unit circle in the plane w±2- 
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We define further the class of compensating factors of 



the simple form [11 



Si m (u) 



1 



(() (m) 



where the exponents, written in terms of the powers 7; 
defined in (j33]) - ([34f and the Kronecker symbol Si m , as 



(44) 



7i°=7i(l + Su), 7^=72(1 + ^2), 



are chosen such that Si m (u) cancel the dominant singu- 
larities defined in (j3"3")l . Following [111 ], we further expand 
the product S l ; m (M)-BRGs(w, y) in powers of the variable 

Wlm(u), as 

Sl m (u)B KGS (u,y) = ^2 C n,RGS 

{y){w lm {u)) n . (45) 

Tl>0 

For the optimal mapping w±2 this expansion converges 
in the whole disk \w±2\ < 1, i.e. in the whole u-plane 
cut along u > 2 and u < — 1, and has the best asymp- 
totic convergence rate [ll|, H3|. Moreover, since the first 
singularities of the Borel transform are compensated by 
the softening factor, a good convergence is expected also 
at finite orders. For other mappings, with either I > 1 
or m > 2, the expansions would converge in the unit 
disks \wi m \ < 1 if the singularities situated between 
—l<u<—l and 2 < u < m were completely removed by 
the compensating factors. In practice, however, only the 
first singularities at u = — 1 and u = 2 are compensated 
by the factors Si m (u), and after the compensation they 
may survive as "mild" branch points, where the function 
vanishes instead of becoming infinite. The presence of 
the residual cuts sets a limit on the convergence domain 
of the expansions (|4"5j) , but for a mild singularity the ef- 
fect is expected to become important only at large orders 

While the optimal mapping is based on a mathematical 
theorem [ll|, H3] , there is no such rigorous result for the 
form of the softening factors. They are arbitrary to a 
large extent and can be chosen empirically. With the 
choice (|43|) . the compensating factors used in different 
expansions are parametrized in different way. By this we 
reduce the bias related to the choice of these factors. 

By combining the expansion (|45[) with the definition 
(|32l) , we are led to the class of BRGS expansions 

Sbrgs(s) = £ cgksd,) W£ GS ( S ), (46) 

n>0 



where 



<RGs(s) = 7T PV / ex P 



U \ {Wl m (u)) r 



AT J \Poa s (~s)J Si m (u) 





and the coefficients c^ r ^ GS (?/) are defined by the expan- 
sion (|4"5]) . 



du, 



(47) 



For completeness we write below the similar "Borel 
and contour-improved" (BCI) expansions @, [ll| 



5bci( S ) = E<c!<ci ) W, 



(48) 



where the expansion functions are expressed in terms of 
the running coupling a s (s): 

oo 

Wn.dW-^Vye S[m{u) du, (49) 



and the coefficients c^™| are defined by the expansion 
S lm (u)B(u) = J2&(^M) n . (50) 



n>0 



Finally, the "Borel improved fixed-order" (BFO) ex- 
pansions are written as p, [ll[ 



5 BF 0( S )=E C n,F0( S )< 



FO' 



(51) 



where the the expansion functions involve the fixed-scale 
coupling a s (M^): 



= ^pv 



and the coefficients c^™ Q (s) are defined by: 

Si m (u)B FO (u, s) = c^fo(s) {wi m {u)) n . (53) 



n>0 



The properties of the new expansions were discussed 
in detail in [H, [2|| in the particular case of the optimal 
mapping w\2 ■ Their definition is based on a prescription 
for the infrared ambiguity of perturbation theory which 
is consistent with analyticity in the energy plane. When 
reexpanded in powers of a s , the expansions reproduce 
the coefficients c„ i known from Feynman diagrams, up 
to any order TV. The remarkable feature is that the ex- 
pansion functions Wn resemble the expanded function 
D(s), being singular at a s = and admitting divergent 
expansions in powers of the coupling. Therefore, the di- 
vergent pattern of the expansion of D(s) in terms of these 
new functions is expected to be tamed. The numerical 
studies reported in the next section confirm this expec- 



tation. As in [ll[, we shall use in these analyses the 



expansions based on the optimal mappings w\2 and the 
alternative mappings w\^ 1 wioo and u>23- 



3 A formal proof for the optimal mapping w\2 was given in |2S( , 
where it was shown that under some special conditions the ex- 
pansion H48H is convergent in a domain of the complex s-plane. 
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VI. MODELS FOR THE ADLER FUNCTIONS 

In order to test numerically the convergence properties 
of the BRGS expansions defined in the previous section, 
we consider a class of physical models discussed recently 
in the literature @, 1, M, M, HI • 

We first consider the model proposed in [BJ , where the 
Adler function Dbj(s) is defined as the PV- regulated 
Laplace-Borel integral 



Po 



exp 



" )B BJ (u)du, (54) 
p a s {-s)J 



in terms of a Borcl transform Bbj{u) parametrized in 
terms of UV and IR renormalons and a regular, polyno- 
mial part: 



B B j{u) _ dUV . 



B^(u) + B^\u) + B^(u) + d™ + d PO ~ 



In the expressions of the renormalons 



(55) 



d m 



l + h(p-u) + 



(p — m)> 
d uv r 



(56) 



most of the parameters were fixed by imposing 
renormalization-group invariance at four loops. Finally, 
the free parameters of the model, namely the residues 
di jdj? and d\ R of the first renormalons and the coef- 
ficients dQ°,di° of the polynomial in (|55p . were deter- 
mined by the requirement of reproducing the perturba- 
tive coefficients c n .i for n < 4 from and the estimate 
c 5: i = 283. Their numerical values are @: 

cfUv = -1.56 x 10 -2 , df = 3.16, 4* = - 13 - 5 > ( 57 ) 

4° = 0.781, 4° = 7.66 x 1(T 3 . (58) 

Then all the higher-order coefficients c„ ; i are fixed and 
exhibit a factorial increase. The numerical values up to 
n = 18 are listed in H, @. 

The magnitude of the residue d 2 R of the first IR renor- 
malon in the above model was questioned by some au- 
thors @, HI|. In order to avoid any bias, we have also in- 
vestigated alternative models, in which a smaller residue 
at u = 2 was imposed (an extreme case of this type of 
alternative models, in which the singularity at u = 2 is 
completely removed, was investigated recently in [32|). 
In one example, we have retained the same expressions 
as in [B[ for the first three singularities and the same val- 
ues of the residues at u = — 1 and u — 3, while choosing 
a smaller residue at u = 2, 4 = 1. In order to repro- 
duce the first five coefficients c„ ; i, the model must con- 
tain then three additional free parameters, which were 
introduced by a quadratic term in the polynomial and 
two additional IR singularities, at u — 4 and u = 5. For 



convenience, the nature of these additional singularities, 
which is not known, was taken to be the same as that 
of the u = 3 singularity. Thus, we have considered the 
alternative model: 

= B^(u) + B^(u) + B^(u)+ d ™ 



(4 - u) 



3.37 



4 R 



iPO 



+ di°u + d 



(5-u) 



3.37 



PO„,2 



(59) 



where, as discussed above, we have taken as input 
4 V ,4 R and 4 R from (|57[) and determined the remain- 
ing five parameters by matching the coefficients c Ui \ for 
n < 5, which gives: 



4 U = 3.2461, d{ 



po 



1.3680, 4° = 0.2785, 



560.614, 



-1985.73. 



(60) 



The physical plausibility of this type of models was 
discussed in several recent works [f| H HH, [32| . In par- 
ticular, arguments in favor of the first model presented 
above were brought in 0, [32| . In the present work we 
adopted these models as a mathematical framework for 
testing the convergence properties of the various expan- 



0.14 




FIG. 1: Real part of the Adler function of the model [f| de- 
fined in (|55[) - H58[) . calculated along the circle s = exp(i#) 
for a s (M?) = 0.3156, using the BRGS perturbative expan- 
sions (|46|l - (|47| l with N = 5 terms. The exact function is 
represented by the solid line. 



We illustrate first the properties of the new expansions 
by the approximation they provide to the expanded func- 
tion along the circle s = M % exp(i#). In Fig. [TJ we show 
the real part of the Adler function for the model Q de- 
fined in (|55[) - (|5"8"|) . calculated with the new BRGS expan- 
sions ()46[) - (l47l) . with N = 5 terms in the perturbative se- 
ries. As in llj, we considered the expansion functions for 
the optimal mapping W12 and the alternative mappings 
W13, u>ioo and W23- For the comparison with previous 
studies @, BEL we have used a s (M^) = 0.3156 in this 
calculation. In Fig(^l we show the same curves for the 
imaginary part of the Adler function. In Figs. [3] and |4j 
we repeat our calculations with N = 18 terms. 
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FIG. 2: As in Fig. [T] for the imaginary part of the Adler 
function. 



0.12 



0.08 



0.04- 




FIG. 3: As in Fig. [T]for N = 18 terms in the expansions. 



As shown in Figs. QJH the BRGS expansions provide 
a good description of the exact function along the whole 
circle, including the points near the timelike axis, which 
correspond to 9 = 0, and near the spacelike axis, where 
8 = ir. The worse approximation provided by the map- 
ping u>23 for N = 18 can be explained by the effect of the 
residual mild cut between u = — 1 and u = — 2, which 
limits the convergence radius of the expansion (|45[) in 



0.08 



0.04 



-0.04 



-0.08 





1.5 2 
9 (rad) 



2.5 



powers of W23 to u < 1. For other mappings, the diver- 
gence due to the residual cuts is manifest only for u > 2, 
and this region is more suppressed by the exponent in the 
Laplace-Borel integrals (|4"T)) defining the expansion func- 



11[). In all the cases, 



FIG. 4: As in Fig. [5] for N = 18 terms in the expansions. 



tions (for more explanations see 
however, the conformal mappings improve the conver- 
gence rate for small values of it, which bring the most 
important contribution to the Laplace-Borel integral. 

Similar representations along the circle \s\ = M% using 
the standard CI and FO expansions are given in Refs. [5|, 
Q . For the standard RGS expansion the results reported 
in [l4| show that its predictions are quite close to those 
of the standard CIPT. The standard expansions exhibit 
bigger and bigger oscillations with increasing N, failing 
to reproduce with accuracy the exact function along the 
circle. 

The BRGS expansions should be compared with the 
Borel improved CI and FO expansions, defined in (|4"5)l - 
(l4"9l a nd (|?T|) - ([5^1) respectively, which were investigated 
in [ll| (the particular expansion based on the optimal 
mapping W12 was treated also in Q ) . As illustrated in the 
figures given in 0, [ll[ , the Borel improved CI expansions 
reproduce very well the exact function, much like the new 
BRGS expansions, while the Borel improved FO expan- 
sions provide a very good approximation near the space- 
like axis, where the powers of the logarithm ln(— s/M^) 
present in the coefficients are small, but the approxima- 
tion becomes worse near the timelike axis, where the log- 
arithm acquires a large imaginary part. 

In order to assess the physical relevance of the con- 
vergence acceleration of the perturbative expansions, we 
considered the behaviour of the new BRGS schemes in 
the context of r-lepton hadronic width, which requires 
the theoretical calculation of the quantity <j'°) defined in 

s(o) 



([3]) . In Tables [Il lIVI we give the differences 
order by order in perturbation theory for the models 
discussed above, using various perturbative expansions. 
The tendency of this quantity to flatten out to would 
indicate that a particular scheme is efficient and reliable. 

In TableJT]we show these differences for the model pro- 
posed in Q and reviewed in Eqs. (|55 ]) -(|55 ]) . and in Table 
|TT]we present the results for the alternative model speci- 
fied in Eqs. ([5i?)) - (|60|) . For a consistent comparison with 
previous results reported in [5], [|| [Til ] , we performed the 
calculations with a s (M%) = 0.34. 

The first three columns of Tables U and |TT] show that 
at low truncation orders the standard FO expansion pro- 
vides a more precise approximation for the model pre- 
sented in Table [U while the standard CI expansion de- 
scribes better alternative models of the type shown in 
Table [TT1 characterized by a smaller residue of the first 
IR renormalon. These features were discussed also in 
@,[33. At lar ger orders, the standard FO expansion ex- 
hibits in both cases a milder divergence, explained [To| 
by the cancellations between the contributions of the co- 
efficients c„.i and the remaining terms in the series (|10|) . 

As concerns the RGS expansion, it provides at low or- 
ders an approximation comparable to the standard CI 
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TABLE I: The difference S (0) -<s££ct for the model B B .j proposed in || and specified in (|55 j) -([58 )l . calculated for a s (M?) = 0.34 
with the standard CI, FO and RGS expansions, and the new BRGS expansions (|46|l - (|47| l for various conformal mappings wi m , 
truncated at order N. Exact value 8^l\ rt = 0.2371. 



N 


CI 


FO 


RGS 


BRGS tow 


BRGS wis 


BRGS toioo 


BRGS w 23 


2 


-0.0595 


-0.0679 


-0.0574 


-0.0347 


-0.0239 


-0.0417 


-0.0177 


3 


-0.0473 


-0.0345 


-0.0440 


-0.0333 


-0.0301 


-0.0349 


-0.0303 


4 


-0.0388 


-0.0171 


-0.0347 


-0.0089 


-0.0142 


-0.0067 


-0.0132 


5 


-0.0349 


-0.0083 


-0.0315 


-0.0070 


-0.0086 


-0.0058 


-0.0070 


6 


-0.0325 


-0.0043 


-0.0284 


-0.0073 


-0.0071 


-0.0064 


-0.0072 


7 


-0.0325 


-0.0029 


-0.0298 


-0.0059 


-0.0057 


-0.0056 


-0.0044 


8 


-0.0354 


-0.0018 


-0.0309 


-0.0041 


-0.0035 


-0.0041 


-0.0011 


9 


-0.0367 


-0.0004 


-0.0363 


-0.0023 


-0.0019 


-0.0028 


-0.0010 


10 


-0.0529 


0.0019 


-0.0483 


0.0014 


-0.0012 


-0.0020 


0.0004 


11 


-0.0409 


0.0031 


-0.0458 


0.0036 


-0.0008 


-0.0016 


-0.0009 


12 


-0.1248 


0.0065 


-0.1335 


0.0031 


-0.0006 


-0.0015 


0.0005 


13 


0.0258 


0.0037 


0.0534 


0.0026 


-0.0004 


-0.0015 


-0.0005 


14 


-0.5286 


0.0204 


-0.7850 


0.0018 


-0.0003 


-0.0015 


-0.0011 


15 


0.8640 


-0.0201 


1.7734 


0.0006 


-0.0002 


-0.0015 


0.0044 


16 


-3.5991 


0.1447 


-7.7043 


0.0001 


-7- 10" 6 


-0.0015 


-0.0131 


17 


9.3560 


-0.4252 


24.8586 


-0.0004 


4 ■ 10~ 6 


-0.0014 


0.0238 


18 


-31.76 


1.907 


-94.26 


-0.0013 


-0.0001 


-0.0013 


-0.0310 



TABLE II: As in Table |l] for the modified model B aH specified in ([59 )l -lf60 |l . Exact value <5^ ct = 0.2102. 



N 


CI 


FO 


RGS 


BRGS W12 


BRGS wi 3 


BRGS wioo 


BRGS w 23 


2 


-0.0326 


-0.0410 


-0.0305 


-0.0078 


0.0030 


-0.0148 


0.0092 


3 


-0.0204 


-0.0076 


-0.0171 


-0.0064 


-0.0033 


-0.0080 


-0.0034 


4 


-0.0119 


0.0098 


-0.0078 


0.0180 


0.0127 


0.0202 


0.0137 


5 


-0.0080 


0.0186 


-0.0046 


0.0199 


0.0183 


0.0211 


0.0110 


6 


-0.0061 


0.0216 


-0.0026 


0.0175 


0.0175 


0.0182 


0.0197 


7 


-0.0061 


0.0188 


-0.0047 


0.0193 


0.0150 


0.0153 


0.0225 


8 


-0.0079 


0.0111 


-0.0052 


0.0201 


0.0132 


0.0131 


0.0258 


9 


-0.0065 


0.0008 


-0.0078 


0.0106 


0.0101 


0.0109 


0.0259 


10 


-0.0178 


-0.0070 


-0.0142 


-0.0012 


0.0047 


0.0083 


0.0273 


11 


0.0022 


-0.0098 


-0.0020 


-0.0118 


-0.0010 


0.0058 


0.0260 


12 


-0.0690 


-0.0031 


-0.0737 


-0.0231 


-0.0054 


0.0037 


0.0274 


13 


0.1019 


0.0015 


0.1397 


-0.0310 


-0.0081 


0.0023 


0.0264 


14 


-0.4207 


0.0242 


-0.6549 


-0.0339 


-0.0093 


0.0014 


0.0258 


15 


1.0234 


-0.0168 


1.9784 


-0.0347 


-0.0086 


0.0008 


0.0313 


16 


-3.3572 


0.1398 


-7.3731 


-0.0316 


-0.0062 


0.0004 


0.0139 


17 


9.7378 


-0.4435 


25.4225 


-0.0239 


-0.0028 


0.0002 


0.0507 


18 


-31.15 


1.874 


-93.316 


-0.0156 


0.0003 


-2 ■ 10~ 5 


-0.0041 



expansion for the first model, and slightly better for the 
second model. However, the description deteriorates be- 
yond N — 10 where large oscillations in the results ap- 
pear, the RGS expansion exhibiting in a more dramatic 
way than CIPT and FOPT the divergent pattern of the 
QCD perturbation theory. An improvement of its large- 



order behaviour by the techniques discussed in this paper 
is therefore mandatory. 

The last four columns of Table U show that for the first 
model the new BRGS expansions provide a very good 
approximation already at low orders, and the accuracy 
increases with the truncation order N. According to the 
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TABLE III: The difference <5 (0) -<j£2«t for the model B B j proposed in @| and specified in l(55 l> -<[58 l> . calculated for a s (M?) = 0.34 
with the improved BCI expansions (1481) - (I49|) and the BRGS expansions (|46[) - (|47| ). for various conformal mappings wi m , truncated 
at order N. Exact value Sl°l rt = 0.2371. 



N 


BCI W12 


BRGS wi2 


BCI wis 


BRGS wi 3 


BCI TOloo 


BRGS toioo 


BCI W23 


BRGS W23 


2 


-0.0394 


-0.0347 


-0.0301 


-0.0239 


-0.0488 


-0.0417 


-0.0248 


-0.0177 


3 


-0.0362 


-0.0333 


-0.0341 


-0.0301 


-0.0396 


-0.0349 


-0.0343 


-0.0303 


4 


-0.0108 


-0.0089 


-0.0177 


-0.0142 


-0.0083 


-0.0067 


-0.0165 


-0.0132 


5 


-0.0081 


-0.0070 


-0.0103 


-0.0086 


-0.0061 


-0.0058 


-0.0079 


-0.0070 


6 


-0.0047 


-0.0073 


-0.0065 


-0.0071 


-0.005 


-0.0064 


-0.0052 


-0.0072 


7 


-0.0032 


-0.0059 


-0.004 


-0.0057 


-0.0038 


-0.0056 


-0.0026 


-0.0044 


8 


-0.0032 


-0.0041 


-0.0028 


-0.0035 


-0.003 


-0.0041 


-0.0024 


-0.0011 


9 


-0.0030 


-0.0023 


-0.0023 


-0.0019 


-0.0025 


-0.0028 


-0.0024 


-0.0010 


10 


-0.0020 


0.0014 


-0.0023 


-0.0012 


-0.0023 


-0.0020 


-0.0018 


0.0004 


11 


-0.0012 


0.0036 


-0.0023 


-0.0008 


-0.0022 


-0.0016 


-0.0023 


-0.0009 


12 


-0.0009 


0.0031 


-0.002 


-0.0006 


-0.0022 


-0.0015 


0.0003 


0.0005 


13 


-0.0009 


0.0026 


-0.0016 


-0.0004 


-0.0022 


-0.0015 


-0.0023 


-0.0005 




-U.UUU / 


0.0018 


-0.001 


-0.0003 


-0.0022 


-0.0015 


0.0024 


-U.UU11 


15 


-0 0004 


0.0006 


-0.0005 


-0.0002 


-0.0021 


-0.0015 


-0.0015 


0044 


16 


-0 0003 


0.0001 


-0.0002 


-7- 10~ 6 


-0.002 


-0.0015 


-0.0028 


-0.0131 


17 


-0 0003 


-0.0004 


0.0001 


4 ■ 10~ 6 


-0.0019 


-0.0014 


0.0162 


0238 


18 


-0 0003 


-0.0013 


0.0002 


-0.0001 


-0.0017 


-0.0013 


-0.0445 


-0.0310 




TABLE IV: 


As in Table ITTT1 for the modified model B a \ t specified in ( 


59j)-d6nj ■ Exact value <5™ c t = 


0.2102. 


N 


BCI wi 2 


BRGS W12 


BCI w 12 


BRGS wis 


BCI w loo 


BRGS w lac 


BCI W23 


BRGS W23 


2 


-0.0125 


-0.0078 


-0.0032 


0.0030 


-0.0219 


-0.0148 


0.0021 


0.0092 


3 


-0.0093 


-0.0064 


-0.0072 


-0.0033 


-0.0127 


-0.0080 


-0.0074 


-0.0034 


4 


0.0161 


0.0180 


0.0092 


0.0127 


0.0186 


0.0202 


0.0104 


0.0137 


5 


0.0188 


0.0199 


0.0166 


0.0183 


0.0208 


0.0211 


0.0190 


0.0110 


6 


0.0161 


0.0175 


0.0169 


0.0175 


0.0182 


0.0182 


0.0158 


0.0197 


7 


0.0099 


0.0193 


0.0118 


0.0150 


0.0128 


0.0153 


0.0072 


0.0225 


8 


0.0100 


0.0201 


0.0062 


0.0132 


0.0080 


0.0131 


0.0034 


0.0258 


9 


0.0073 


0.0106 


0.0041 


0.0101 


0.0052 


0.0109 


0.0036 


0.0259 


10 


-0.0047 


-0.0012 


0.0042 


0.0047 


0.0044 


0.0083 


0.0013 


0.0273 


11 


-0.0120 


-0.0118 


0.0034 


-0.0010 


0.0044 


0.0058 


-0.0034 


0.0260 


12 


-0.0095 


-0.0231 


0.0009 


-0.0054 


0.0046 


0.0037 


-0.0021 


0.0274 


13 


-0.0080 


-0.0310 


-0.0016 


-0.0081 


0.0047 


0.0023 


-0.0042 


0.0264 


14 


-0.0101 


-0.0339 


-0.0028 


-0.0093 


0.0044 


0.0014 


0.0022 


0.0258 


15 


-0.0093 


-0.0347 


-0.0023 


-0.0086 


0.0040 


0.0008 


-0.0015 


0.0313 


16 


-0.0058 


-0.0316 


-0.0011 


-0.0062 


0.0034 


0.0004 


-0.0029 


0.0139 


17 


-0.0043 


-0.0239 


0.0 


-0.0028 


0.0028 


0.0002 


0.0173 


0.0507 


18 


-0.0044 


-0.0156 


0.0005 


0.0003 


0.0022 


-2 ■ 10" 5 


-0.0485 


-0.0041 


recent work 32 1, this model is a solid candidate for the 
physical Adler function. For the alternative model, the 
results shown in Table ITTI indicate a slightly worse approx- 
imation at low orders. However, the good convergence of 
the new expansions at large orders, in contrast with the 
big divergencies of the standard expansions, is visible also 


in this case. 

As we mentioned, the standard RGS expansion is 
rather similar to the standard CI expansion up to rel- 
atively large orders, beyond which the RGS expansion 
starts to exhibit much wilder oscillations. It is of interest 
to compare these schemes also in the Borel improved ver- 
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sions given in Eqs. (|3S ]) -([i7 |) and ([I8" |) -(|4"9" ]) . respectively. 
This comparison is presented in Tables IIIII and IIVI The 
results show that, with small variations, the approxima- 
tion provided by the BCI and BRGS expansions is similar 
up to large truncation orders N, for both models consid- 
ered. Thus, using the technique of series acceleration by 
conformal mappings of the Borel plane, the strong diver- 
gence of the standard RGS was considerably tamed. We 
recall that an advantage of the BRGS expansion is that 
it does not require the numerical determination of the 
running coupling along the integration circle \s\ = M 2 , 
involving only analytical expressions. 

VII. DETERMINATION OF a s (M?) 

In this section we shall use the new BRGS expansions 
defined in the present paper for a new determination of 
a s (M 2 ) in the MS scheme. The determination of this 
fundamental parameter is one of the important goals 
of this work and significant care has to be exercised in 
adopting proper values of input along with the experi- 
mental and theoretical uncertainties. 

We use as input the recent phenomenological value of 
the pure perturbative correction to the hadronic r width 

m 

6 SL = °- 2037 ± 0.0040 cxp ± 0.0037pc, (61) 

where the first error is experimental and the second re- 
flects the uncertainty of the higher-order terms ("power 
corrections") in the OPE estimated by reasonable theo- 
retical assumptions. We emphasize that our calculation 
is not based on the models discussed in the previous sec- 
tion, but relies only on the known coefficients c n ^\ given 
in and the conservative choice Cs^ = 283 ± 283 for 
the next coefficient [!, [l(| ■ 

Using this input, the values of a s (M 2 ) obtained with 
the new BRGS expansions defined in (l46l) - (l47l) . with 

the expansion functions W^ GS , W^ GS , and 
RGS ' respectively, arc: 

0.3189 ± 0.0034 cxp ± 0.0031p C tomtlM to.^isO^), 
0.3198 ± 0.0034 cxp ± 0.0031 PC lo.rass^i) -omoUM, 
0.3180 ± 0.0034 cxp ± 0.0031p C tffitei) ±SS(&), 
0.3188 ± 0.0034 oxp ± 0.0031 PC jftgl$(c 6 i) l&oooS^). 

(62) 

The first two errors are due to the uncertainties of the 
phenomenological value of given in (pjTj) . The third 
one, produced by the range adopted for the coefficient 
C5 ; i, brings the most important contribution to the to- 
tal error. The last uncertainty accounts for the higher 
terms in the expansion of the /3 function, simulated as 
in by an additional coefficient ,84 = ±/3f//?2 in this 
expansion. We explored also the influence of the scale 



variation, choosing it as fi 2 = £M 2 with £ = 1 ± 0.5, 
but the effects are very small, so we show only the result 
corresponding to the choice = M 2 in ((26)) . 

A very small sensitivity of a s (M 2 ) to the variation of 
the scale is specific also to the standard CIPT analyses 
0[§], the RGS expansion (l4| . and the Borel improved 
CI expansions @, [Tl[ . The uncertainty related to the co- 
efficient C54 is bigger in the case of the Borel improved 
expansions than in the standard CI and RGS. However, 
as discussed in Jf| [ll| , having in view the divergent char- 
acter of the series, the truncation error in the latter ver- 
sions is certainly underestimated. The calculation of the 
five-loop coefficient 05,1 is therefore of great interest, as it 
would reduce considerably the total error of a s (M 2 ) de- 
termined from the Borel improved perturbation schemes. 

By taking the average of the central values and of the 
errors given in Eq. (|62p we obtain the prediction 

a s (M T 2 ) = 0.3189 ±0.0034 oxp ±0.0031p C (63) 

±o:oioi(c 8 ,i) ±o.oo% 4 , 

which becomes, after adding the errors in quadrature, 

a s (M r 2 ) = 0.3189 ±°;°i« . (64) 

We emphasize that the error quoted above was ob- 
tained as the average of the errors of the individual 
determinations ([52"]) . Much lower uncertainties would 
have been obtained if standard statistical procedures for 
combining independent determinations were applied. In 
practice, although the values given in (|6"2")l may be con- 
sidered independent theoretical determinations, we pre- 
fer the conservative errors given in (|63[) . which avoid 
any bias. We note however the remarkable consistency 
of the theoretical determinations given in (|62p . which 
is a strong argument in favor of our predictions. It is 
remarkable also that the central value of our predic- 
tion (fS4")) practically coincides with the world average 
a s (M T 2 ) = 0.3186 ± 0.0056 Q. 

By evolving (1641) to the scale of Mz, using the CRun- 
Dec package [44j . our prediction reads 

a s {M%) = 0.1184 (65) 

where the central value coincides with the 2012 world 
average, a s (M§) = 0.1184 ± 0.0007 Q. 

It is of interest to compare the result (|6~4")l obtained 
with the BRGS expansions defined in the present work 
with other recent determinations of a s (M 2 ). The values 
reported in the recent works 0-[3 are not all based on 
the same input. Therefore, for a consistent comparison of 
different perturbative schemes, we use in what follows the 
same input for the phenomenological value of <5 p ^ en , given 
in (|6"Tj) . and the five-loop coefficient, C54 = 283 ± 283. 
Then the standard FO, CI and RGS expansions lead to 
the predictions 

a s (M 2 ) = 0.3199 FO , 

a s (M 2 ) = 0.3419 ± ;° 085 CI, 

a s {M 2 ) = 0.3378 ±g;°g|| RGS, (66) 
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where the FO result is quoted in [T(| and the RGS one in 
[3] ■ The FO and CI expansions improved by conformal 
mappings of the Borel plane (denoted as BFO and BCI, 
respectively) give [EEI 

a s (Af T 2 ) = 0.3109 BFO, 
a s (M r 2 ) = 0.3195 ±8;8i|» BCI. (67) 

One can see that the BRGS and BCI expansions, which 
implement both Borel and RG-summation, lead to simi- 
lar values for a s (M 2 ), which actually are also close to the 
prediction of standard FOPT, where neither of the two 
summations is performed. On the other hand, the stan- 
dard CI and RGS expansions, where only the summa- 
tion related to RG is implemented, lead to values larger 
by about 0.02, while the BFO expansions, which im- 
prove only the large-order behaviour without summing 
the terms known from the RG, lead to a value smaller 
by about 0.01 than the predictions of BRGS, BCI and 
FOPT. 

In the recent work, Ref. [l3|, it has been pointed out, 
on the basis of a detailed analysis of moments of the spec- 
tral functions measured by OPAL experiment, that the 
uncertainty of the nonperturbative contributions to the 
hadronic T-width may be larger than usually assumed. 
To account for this possibility, we adopt a more conser- 
vative value for the uncertainty of S^ cn due to power 
corrections, viz. replacing ±0.0037 in (|6"Tj) by the esti- 
mate ±0.012 of this uncertainty quoted in (l3| . 

Using this input, the value ±0.0031 of the error on 
a s (M^) produced by the power corrections, quoted in 
(|63| . increases to i" 'oi03 > leading to 

a s {M 2 T )= 0.3189 ±°;°15?- (68) 

We shall adopt this result, having the same central value 
and a more conservative error than in (|64p . as our final 
prediction. 

By evolving (|55|) to the scale of Mz our prediction 
reads 

a s (M|) = 0.1184 ±° Zl (69) 

VIII. DISCUSSION AND CONCLUSIONS 

The non-strange hadronic decays of the r lepton pro- 
vide one of the most important ways of extracting the 
strong coupling a s . The pcrturbative schemes of the 
Adler function that are used in this extraction continue, 
however, to be a significant source of uncertainty. There 
are two major ambiguities that affect the QCD perturba- 
tion expansions: one is related to the implementation of 
the renormalization-group invariance, the other regards 
the large-order behaviour of the series. 

The first ambiguity is usually illustrated by the signifi- 
cant difference between the predictions of the FO and CI 
expansions of the Adler function. In a recent work, Ref. 



[14} . we used for the analysis of the r hadronic width 
another renormalization-group-improvcd expansion, pro- 
posed in [13, Gil ! which sums in an analytically closed 
form the logarithms accessible by RG invariance. 

The second ambiguity is associated with the fact that 
the perturbative series are divergent, due to a factorial 
growth of the coefficients calculated from Feynman di- 
agrams [H-[22j. In QCD this ambiguity is even more 
dramatic than in QED due to the IR renormalons, sin- 
gularities in the Borel plane situated on the positive axis, 
which prevent the unique reconstruction of the function 
by means of a Laplace-Borel integral. However, if one 
adopts a suitable prescription (principal value), it is pos- 
sible to exploit the available knowledge on the large-order 
behaviour of the coefficients for defining a new expansion, 
in which the divergent pattern is tamed to a great extent. 
Such a method was proposed some time ago in [27)429] 
and was applied recently in 0, [||, [ll| to both the FO and 
CI expansions of the Adler function. It is based on series 
acceleration by analytic continuation in the Borel plane, 
achieved by conformal mappings and the "softening" of 
the dominant singularities of the Borel transform. 

In the present work, we used this method in order to 
improve the large-order behaviour of the RGS expansion 
discussed in [l4j. That it should be possible to have a 
straightforward application of the technique is not obvi- 
ous, as the RGS expansion ([T2"]) involves a set of com- 
plicated functions D n (y) determined iteratively by the 
differential equations (|15[) . However, after expressing the 
series in the alternative form (|25[) .ie.. as an expansion in 
powers of the one- loop coupling (|2l))) , it was possible to 
show in Sec. IIVI that the dominant singularities of this 
expansion in the Borel plane coincide with those of the 
standard Borel transform B(u) of the CI expansion. As 
a result, we are able to apply the techniques discussed 
in [H, [H, [ll| , defining the class of improved expansions 
(gnjl-dlZl), denoted as BRGS expansions. 

As discussed in earlier works [ill, HH , the expansions 
based on conformal mappings of the Borel plane have a 
number of remarkable properties. In particular, the di- 
vergent pattern of the expansion of D(s) in terms of these 
new functions is expected to be tamed. The detailed nu- 
merical studies of two representative models for the Adler 
function, presented in Sec. IVI[ show that indeed the ex- 
pansions improved by both RG summation and analytic 
continuation in the Borel plane, i.e. BCI and BRGS, 
provide a very good approximation of the true functions 
in the complex energy plane up to high orders. This is 
seen from Figs. [THof Sec. IVII and the similar figures 
presented in Refs. [fl [TTJ] . As a consequence, as shown 
in Tables U - IIVI the approximation of the integral 5^ 
provided by these expansions is very good and improves 
with increasing N. In contrast, the standard expansions 
shown in the first three columns of Tables HI and ITTI exhibit 
a divergent pattern. 

As argued in [f| [l(| and is seen also from our results, 
the FO expansion of the observable 5^ exhibits a better 
behaviour at high orders, which results from suitable can- 
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collations between the increasing coefficients c n ^\ and the 
remaining terms in (|10|) . This cancellation is spoiled in 
the standard CI expansion (fTTj) . which sums only the RG 
part, leaving the increase of c n> \ uncompensated. The 
better large-order behaviour of FOPT is one of the argu- 
ments used in favor of this expansion compared to CIPT 
0. As concerns the standard RGS expansion, some 
cancellations are taking place at low orders, but above 
N = 10 this expansion exhibits a divergence pattern even 
more dramatic than CIPT. 

It should be noted however that the better behaviour 
of FOPT is restricted only to integrals like 5^ or some 
special moments of the spectral function: as shown in 
previous studies of specific models @, [ll[ , the point- 
wise description of the true Adler function along the cir- 
cle \s\ = M% provided by the FO expansion is quite poor. 
The good convergence obtained for some integrals along 
the circle is due to fortuitous cancellations of large contri- 
butions from different integration regions. On the other 
hand, as we mentioned, the BCI and BRGS expansions 
yield a good pointwise approximation of the true function 
along the whole circle. 

For the determination of a s (M^) presented in Sec. IVIII 
we used as input the phenomenological value of the QCD 
correction 6^ to the t hadronic width evaluated in [l(| . 
Our analysis shows that the BRGS and BCI expansions, 
improved by both RG and high-order Borel summation, 
lead to similar results for a s (M^), which are actually 
also very close to the standard FOPT prediction obtained 



with the same input. It is remarkable that the central 
value of the prediction (|M|) obtained with the new BRGS 
expansions coincides with the world average quoted in [l[ . 
We also considered the possibility that the uncertainty 
related to nonperturbativc contributions might be larger, 
as follows from the recent analysis [l3[ . The error of our 
final prediction (|68p is a conservative estimate that takes 
into account this possibility. 

In conclusion, we advocate the use of the BCI and 
BRGS expansions of the QCD Adler functions, which im- 
plement simultaneously the RG invariance and the avail- 
able knowledge about the large-order behaviour of per- 
turbation theory. In particular, the BRGS expansions 
proposed in the present work have the advantage that RG 
summation is implemented through analytically closed 
expressions. Therefore, these expansions are suitable for 
more detailed analyses of r hadronic decays, based on 
the moments of the hadronic spectral function. 
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